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QCD Down Under: Building Bridges 

M.R. Pennington^ 

'^Institute for Particle Physics Phenomenology, University of Durham, Durham DHl 3LE, U.K. 

The strong coupling regime of QCD is responsible for 99% of hadronic phenomena. Though considerable 
progress has been made in solving QCD in this non-perturbative region, we nevertheless have to rely on a 
disparate range of models and approximations. If we are to gain an understanding of the underlying physics and 
not just have numerical answers from computing "black" boxes, we must build bridges between the parameter 
space where models and approximations are valid to the regime describing experiment, and between the different 
modellings of strong dynamics. We describe here how the Schwinger-Dyson/Bethe-Salpeter approach provides 
just such a bridge, linking physics, the lattice and experiment. 



1. LIGHT QUARK MASSES 

This meeting which starts today in the beauti- 
ful Barossa valley is about "building bridges" be- 
tween different approaches to QCD, and between 
these theoretical studies and experiment. I will 
leave others to discuss work on the light cone and 
on flux tubes. Consequently, this talk will serve 
as an introduction to that by Peter Tandy 0, 
which follows. He will provide the details. I will 
give the broad brush description and introduce 
the basic ideas. 

Most discussions of QCD begin with the per- 
turbative calculations which are so successful in 
describing a whole range of hadronic phenomena 
from deep inelastic scattering to beauty produc- 
tion at the Tevatron. These calculations work be- 
cause they study very short distance interactions 
well inside the femto-universe when the vacuum 
appears to be essentially empty, like that of QED. 
However, what makes QCD so much more fasci- 
nating is the fact that over the distance of a fermi 
it becomes strong. Strong physics is responsi- 
ble for confinement, for chiral symmetry breaking 
and for the whole spectrum of hadrons. This is 
what I want to review in this talk. The vacuum is 
then far from empty. It is not just a sea of qq pairs 
and a cloud of gluons, the effects of which can be 
perturbatively computed, but so strong is the in- 
teraction that the quarks, antiquarks and gluons 
form condensates that change the very nature of 



the vacuum. The effect of this is illustrated by the 
way the mass of an up or down quark changes as it 
propagates over the size of a hadron. As pictured 
in Fig. 1, at short distances the u/d quark has a 
very small current mass but over bigger distances, 
the size of a fermi, the medium through which 
it travels generates a constituent mass, which is 
about a third of the mass of a proton. 




r(fm) 

Figure 1. Schematic illustration of how the mass 
of an up and down quark depends on the distance 
r over which it propagates. The quark becomes 
increasingly dressed from 'current to constituent' 
as r increases. 
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Figure 2. Schwinger-Dyson equations for the 2 
and 3-point functions of QED. The quantities 
marked with a sohd dot are dressed. The sohd 
Unes represent fermions, the wiggly hues photons. 



Like aU strong physics problems, this requires 
calculations beyond the perturbative. 

Strong physics problems dominate the world 
of light quarks, (ups and downs). The fact that 
these quarks are light is crucial to the real world. 
Indeed hadron jets in e' annihilation remem- 
ber the direction and spin of the initially pro- 
duced quark and antiquark at a centre-of-mass 
energy as low as a few GeV entirely because most 
hadronisation requires the creation of pairs 
with current masses of just a few MeV. These 
in turn create pions, which are so much lighter 
than any other hadron. If the lightest quark were 
1.5 GeV in mass, like the charm quark, the world 
would be quite different and we would have lit- 
tle hint in hadron interactions of the underlying 
quark dynamics till very much higher energies. 

But how does this picture, shown in Fig. 1, 
emerge from QCD. To answer this, let us consider 
more generally when can masses be dynamically 
generated in a gauge theory — can the bare mass 
be zero and yet the particle have mass [2E1? It is 
well known that in perturbation theory, the mass 
of any particle is proportional to its bare mass 
at every order. So if the bare mass is zero, the 
dressed mass remains zero. Consequently, mass 
generation must be a strong physics problem. A 
possible approach to studying this would be to 
consider the field theory on the lattice. However, 



massless particles do not fit on a finite size lat- 
tice, so calculations have to be performed with 
non-zero mass and the massless result obtained 
by extrapolation. Such an extrapolation has lit- 
tle to do with the lattice computations, but has 
to be calculated in some other way, as we will 
discuss again later. This means we have to treat 
massless gauge theories in the continuum. The 
field equations of the theory are the Schwinger- 
Dyson equations. Fig. 2. These are genuinely 
non-perturbative and are in the form of a set of 
nested integral equations. 

To study the question of dynamical mass gen- 
eration at its simplest, let us consider the elec- 
tron in QED. The dressed fermion propagator is 
determined from the bare fermion propagator by 
the dressed fermion, the dressed photon and the 
full fermion-photon interaction, as in Fig. 2. The 
dressed photon is in turn determined by dressed 
fermions and the full fermion-photon vertex. One 
can imagine solving the coupled electron and pho- 
ton 2-point equations, if we know the full fermion- 
photon interaction, but this of course satisfies its 
own Schwinger-Dyson equation (Fig. 2) that re- 
lates it to the 2 and 3-point functions we have 
discussed and a new 4-fermion interaction, which 
in turn is related to 2, 3, 4 and 5-point functions. 
We have an infinite system of coupled equations. 
It appears that to find the 2-point function we 
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Figure 3. Schwinger-Dyson equation for the 
inverse fermion propagator in an approxima- 
tion, where only the fermion 2-point function is 
dressed. This corresponds to the summation of 
the rainbow graphs shown in the lower equation. 
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Figure 4. Fermion mass defined at Euclidean mo- 
menta by AA{p^ = m?) = m in units of the ul- 
traviolet cut-off, function of the constant 
coupling ao, with zero bare mass, showing how a 
dynamical mass is generated for > — 7r/3. 

need to know all n-point functions, unless we can 
find a way to truncate this system. The only 
self-consistent truncation we know is perturba- 
tion theory, but this is quite inappropriate for 
the problem of mass generation. So let us in- 
stead make a simple, but brutal truncation of the 
Schwinger-Dyson equations. 

Let us treat the full fermion-photon interac- 
tion as bare, viz. 7'', and quench the photon 
propagator in Fig. 2. This gives the rainbow ap- 
proximation shown in Fig. 3. Then the coupling 
a is a constant, ao, that does not run. Conse- 
quently, the only momentum scale is provided by 
the ultraviolet cut-off k, which would be replaced 
by some physical scale on renormalization. With 
this butchery we have a non-perturbative equa- 
tion for the electron propagator, which couples its 
wavefunction renormalization and its mass func- 
tion. In the Landau gauge the wavefunction is 
unrenormalised and we just have to solve for the 
electron mass. An infinite set of equations has 
been reduced to one. This can then be readily 
solved |2I3| . One finds, as shown in Fig. 4, that 
with a zero bare mass the dressed electron re- 
mains massless until the coupling becomes strong 
at some critical value, which in this case is 7r/3. 



As the coupling strengthens beyond this value, a 
larger and larger mass is generated. What this 
means in the world of quenched QED is illus- 
trated by considering the behaviour of a mass- 
less electron in the presence of nuclei of increas- 
ing atomic number. The electron would move at 
the speed of light until the atomic number was 
greater than 140 and then it would slow up as 
the strength of the interaction generates a mass. 

If one performs this calculation in any other 
covariant gauge than Landau a mass is gener- 
ated above some critical value, but this coupling 
strength depends on the gauge, which is of course 
unphysical. This is because the rainbow approxi- 
mation violates the Ward-Green- Takahashi iden- 
tity (WGTI) ^ . To make the answer gauge inde- 
pendent, one not only has to make the interaction 
fulfill the WGTI but the fermion equation must 
be multiplicatively renormalisable too This 
requirement essentially determines the projection 
of the full fermion-photon interaction imposed by 
the fermion Schwinger-Dyson equation — at least 
in quenched QED. We learn that if the coupling 
is strong enough then masses can be generated 
and that is all we need to motivate a programme 
of study in QCD, where the u and d quarks are 
so nearly massless, but hadrons have mass! 




Figure 5. Schwinger-Dyson equations for the 
quark and giuon propagators in QCD in a co- 
variant gauge. The solid lines are fermions, the 
wiggly lines gluons and dashed ghosts. The solid 
dot indicates fully dressed quantities. 
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Figure 6. Coupling in QCD as a function of mo- 
mentum, p. For p 3> ^QCD the running shows 
asymptotic freedom, while for p < Aqcd we 
have different possibilities of confining behaviour. 
Which is the one required by QCD is an issue of 
current research. 




p' (GeV') 

Figure 7. Euclidean mass function Ai{p^) as a 
function of momentum, p, for different current 
masses corresponding to b, s and u, d quarks. The 
chiral limit of massless current quarks is shown by 
the circles, from the calculation of Ref. |12j . 



QCD is inherently more complicated. There 
are of course the gluon self-interactions as well 
as ghosts in covariant gauges. Fig. 5. The first 
Schwinger-Dyson studies of the infrared coupling 
in QCD were by Pagels [H], Mandelstam [Jj, and 
Bar-Gadda "H", ^ with an extensive numerical in- 
vestigation by Nick Brown and myself m. Here 
the ghosts were treated perturbatively to pro- 
duce the right behaviour in the ultraviolet regime. 
The solutions of the equations of Fig. 5 revealed 
a gluon propagator that in the infrared became 
strongly enhanced, as did the strong coupling like 
the upper curve in Fig. 6. This enhancement was 
consistent with a linear confining heavy quark po- 
tential and strong enough to produce a dynamical 
mass for the u and d quarks. This was studied 
in detail by Maris and Roberts who investi- 
gated the quark equation for different bare masses 
from zero to 5 GeV to span the range from up 
and down to beauty, Fig. 7. Provided the com- 
bined effect of the quark-giuon interaction and 
the gluon propagator over the infrared regime is 

^other pioneering studies like that of Baker, Ball and 
Zachariasen |5| were in axial gauges, but they used a trun- 
cation that was subsequently shown to be seriously incon- 
sistent 1101 . 



enhanced, a dressed mass of 350-500 MeV is gen- 
erated, as phenomenology requires — the scale 
being set by Aqcd- The difference between the 
behaviour with a massless current quark and one 
of 3-5 MeV is seen to be very small — except on 
the logarithmic scale of Fig. 7! 

Dynamical mass generation requires the devel- 
opment of non-zero condensates. Their value is 
only strictly defined in the chiral limit. The 
asymptotic behaviour of the mass function is 
given by 

Mij,^) ^ ,n, (l.$)" + C<I|> (1.4)"" 

where d is the appropriate anomalous dimension 
and mo is the current mass. In the limit mg — > 0, 
( qq ) becomes gauge invariant. By looking at the 
ultraviolet behaviour in Fig. 7, Langfeld et al. [T^ 
show a qq condensate of scale —(250 — 300MeV)^ 
results, just as the phenomenology of QCD sum- 
rules |14| and experiments probing low energy tttt 
interactions |15j would claim. 

Further studies have however indicated that 
ghosts might well play a more important role than 
hitherto expected. This follows from a systematic 
investigation by the Tiibingen group of the 
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Figure 8. Momentum dependence of the Landau 
gauge gluon and ghost dressing functions from 
Schwinger-Dyson studies flB^. These give the ef- 
fective quark-gluon coupHng shown. 



coupled ghost, gluon and then quark Schwinger- 
Dyson equations. With no equivalent of the 
Slavnov- Taylor identity to guide the construction 
of the needed ghost-gluon vertex these studies 
have so far been restricted to the Landau gauge, 
where the ghost kernels are known to be sim- 
ple Alkofer and co-workers found that 
contrary to previous investigations the gluon does 
not become enhanced at low momenta, but in fact 
its dressing function goes to zero. In contrast, the 
ghost dressing function is the one that is infrared 
enhanced. The combination of these means that 
the effective coupling does still increase at smaller 
momenta, but with a constant infrared limit, as 
shown in Fig. 8. Though the exact value of this 
limit is known to be dependent on the approxi- 
mations made, the fact that it is larger than 1 is 
sufficient to ensure the dynamical chiral symme- 
try breaking indicated by experiment. 

This scheme appears to solve one mystery: are 
gluons confined or confining 18 ? If a particle 



is confined then one does not expect it to have a 
mass-shell and so the numerator of its propagator 
should cancel any potential pole in the denomi- 
nator. This the Tiibingen scheme J16| achieves 
for the transverse gluons. Nevertheless, we know 
the effective coupling must be enhanced in the in- 
frared if confinement is to be generated. Here this 
is produced by the unphysical degrees of freedom 
of the gluon (and ghost) propagators. Thus this 
scheme for the ghost/gluon sector potentially re- 
solves the dilemma and both confines and is con- 
fined. However, how this mechanism can actually 
generate a near linearly confining potential for 
heavy quark systems is an important issue cur- 
rently under detailed investigation. 

2. LIGHT HADRON MASSES 

One can input the form of the gluon, ghost 
and quark functions just discussed into the bound 
state equations and study the properties of the 
meson spectrum, particularly light hadrons for 
which the long range nature of the forces is so 
critical [Ij. 

As emphasised by Roberts and collaborators 
|19| . the axial Ward identity ensures that the qq 
bound state with pseudoscalar quantum numbers 
is a Goldstone boson with its interactions gov- 
erned by PCAC. In contrast the bound states 
with scalar and vector quantum numbers have 
masses reflecting the mass of the fully dressed 
(or constituent) quark. The scalar mass is found 
to be highly sensitive to the details of the quark 
scattering kernel [201 of Fig. 9, the vector me- 
son states much less so. Consequently, we will 
consider these again shortly. Peter Tandy ^ 
will flesh out many of these considerations when 
he describes his extensive study of bound states. 
The behaviour of the gluon and ghost propaga- 
tors built into these calculations can be com- 
pared with Monte Carlo lattice simulations and 
are in excellent agreement [21116) . While lat- 
tice calculations can only be performed with size- 
able quark masses, the Schwinger-Dyson/Bethe- 
Salpeter system being continuum equations can 
be computed in the massless limit. 

How to continue lattice results to small quark 
masses has been a major research theme over the 
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Figure 9. Bound state equations of Bethe and 
Salpeter for pseudoscalar and vector mesons. 
Their solution depends critically on the qq scat- 
tering kernel indicated in hatched grey. 

last couple of years For instance, the p- 

meson mass depends linearly on the quark mass, 
with a non-zero value in the chiral limit. For large 
mass the pion too depends linearly on the quark 
mass, but its Goldstone nature means it must 
have vanishing mass when its quark components 
are massless. Indeed, if chiral symmetry break- 
ing is dominated by a non-zero gg-condensate 
then the pion mass depends on the square root 
of the quark mass at small values. Lattice results 
with quark mass larger than 150 MeV confirm 
the linear relation between the p and tt masses 
with barely any deviation. In contrast, chiral per- 
turbation theory (^PT) requires the square root 
mass dependence and the presence of logarithms 
of the pion mass. The key to extrapolation is to 
know when these set in and how to match these 
to the behaviour found by lattice calculations. 

If xPT is computed in dimensional regulari- 
sation, this is problematic since infrared and ul- 
traviolet behaviours become entangled. Thus as 
noted by Donoghue et al. these effects in- 
crease at larger pseudoscalar mass, but larger 
masses define the domain of momentum the lat- 
tice should treat exactly. The need is to sepa- 
rate the infrared behaviour where chiral dynam- 
ics dominates and the lattice treats poorly, from 
larger momenta where the lattice embodies the 
correct physics. Thomas et al. |21] from Adelaide 
and Donoghue et al. [23 are amongst several sets 
of authors who have given prescriptions for han- 
dling this transition. 



Here I want to advertise that the Schwinger- 
Dyson/Bethe-Salpeter approach naturally en- 
compasses physics at all momentum scales. Re- 
specting the axial Ward identities ensures that 
chiral dynamics is included with no need to spec- 
ify at which momentum this chiral behaviour 
dominates. The bound state equations determine 
the dependence of the hadron masses for all quark 
masses, which would be exact if we knew the qq 
scattering kernel precisely. In principle lattice 
calculations do the same. In practice they do so 
only for sizeable quark mass. Indeed, the kaon 
mass is the lightest pseudoscalar mass that com- 
putational limitations currently allow. On the 
lattice all masses are only specified in units of lat- 
tice spacing. One can eliminate this dependence 
entirely by plotting the vector meson mass, for 
instance, as a function of the pion mass. Results 
from CP-PACS for different lattice couplings 
are shown in Fig. 10. 

Using the Tiibingen modelling of the strong 
coupling limit of QCD, we can then calculate the 
very same relationship for all quark masses in the 
continuum approach. Following the work of Maris 
and Tandy [27] , the dressed gluon propagator can 
be represented by 

A^-^ip) = C 4 exp ^o'iP) 

where Aq is the bare transverse gluon propaga- 
tor. The parameter C is related to the strength 
of the interaction, while u corresponds to the mo- 
mentum at which the gluon dressing peaks shown 
in Fig. 8 gHI- From the work of Ref. E?] this is 
known to be around 0.5 GeV. Using the simple 
rainbow-ladder approximation to solve the Bethe- 
Salpeter equation for the p and 7r-mesons one 
finds the behaviour shown as the line in Fig. 10. 

Optimal agreement [2^1 between Bethe- 
Salpeter (BS) and CP-PACS results occurs for 
uj — 0.425 GeV, very close to the expected value. 
The continuation from a pseudoscalar mass of 
500 MeV down to a 140 MeV and the chiral limit 
is specified. 

For this continuation one need not estimate 
where the transition to the dominance of chiral 
dynamics occurs. This is built in. This non- 
perturbative approach to the continuum knows. 
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Figure 10. Vector meson mass as a function of 
pseudoscalar meson mass. The data points are 
from the lattice calculations of CP-PACS as 
shown in the inset. The curve that connects the 
crosses is the result of the Bethe-Salpeter calcu- 
lation using the Tiibingen modelling of the full 
quark and gluon functions by Watson et al. [2^1 
The parameter ui is explained in the text. 



However, at the present level of approximation 
this may not yet accord precisely with xPT. Cal- 
culation does show that the pion mass is propor- 
tional to the square root of the current quark 
mass. However, how much of the detailed be- 
haviour embodied in xPT encoded in the chiral 
logs is there in the present solutions of the bound 
state equations is under study [201 • The calcula- 
tions compared in Fig. 10 in the continuum and 
on the lattice are both performed with just two 
quark flavours. The strange quark is ignored. Of 
course, at large pion mass (i.e. greater than 450 
MeV) the p cannot decay. It is stable. The Bethe- 
Salpeter calculation by Watson et al. [2^1 does not 
yet include the contribution from pion loops ei- 
ther. This is known to give a shift of 60 MeV 
or so in the p mass if it is to generate the exper- 
imental width, but elements of double counting 
occur that are under study. Nevertheless progress 
surely has been made. It is not that wc have lat- 



tice calculations valid for larger quark mass and 
chiral perturbation theory at small quark mass 
with educated guesses of how to join them. The 
Schwinger-Dyson/Bethe-Salpeter approach in the 
continuum holds out the prospect that we have 
results applicable at all quark masses. 

In this talk and that that follows by Peter 
Tandy 1 , we see that a bridge is being built 
that not only relates theory to experiment, but 
can relate the lattice to the continuum. The 
Schwinger-Dyson/Bethe-Salpeter connection pro- 
vides a wide-spanning bridge. The other talks in 
the week ahead will doubtless reveal a panorama 
of interlinking bridges between hadron physics 
and the fascinating strong coupling regime of 
QCD down under, that is the subject of this ex- 
citing Workshop. 
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